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ABSTRACT 

We demonstrate that the generalization of the relativistic Toda chain (RTC) is a special reduction of two- 
dimensional Toda Lattice hierarchy (2DTL). This reduction implies that the RTC is gauge equivalent to the 
discrete AKNS hierarchy and, which is the same, to the two-component Volterra hierarchy while its forced (semi- 
infinite) variant is described by the unitary matrix integral. The integrable properties of the RTC hierarchy 
are revealed in different frameworks of: Lax representation, orthogonal polynomial systems, and r-function ap- 
proach. Relativistic Toda molecule hierarchy is also considered, along with the forced RTC. Some applications to 
biorthogonal polynomial systems are discussed. 
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1 Introduction 



Since the paper of Ruijsenaars , where has been proposed, the relativistic Toda chain (RTC) system was investigated in many 
papers [][][] This system can be defined by the equation: 

ci , • V1 i • \ ex P(?"+i - In) 

q n = (1 + eg w )(l + eq n+ i) - 

1 + e^exp(q n+1 - q„) ^ 

,, . ■ w-, , . s exp(<j n - 

- (l + eg„_i)(l + eg„)— — 

1 + e A exp(q n - g„_ij 

which transforms to the ordinary (non-relativistic) Toda chain (TC) in the evident limit e — + 0. The RTC is integrable, which 
was discussed in different frameworks (see, for example, [] [j Q [] an d references therein). In particular, in Q there was 
introduced a Lax triad for the RTC, in the paper |jj there were investigated Lax representation and its orbit interpretation for 
the RTC and the relativistic Toda molecule, and also suggested a discrete version of the RTC which has the same integrals 
of motion (i.e. the same Lax operator, but different evolution). At last, authors of presented the bilinear (Hirota type) 
equations whose solutions satisfy the RTC, implying the possibility of r-function interpretation of the model. Thus, the RTC 
can be obtained from many different starting points. In particular, this is a limit of the general Ruijsenaars system jlj. 

On the other hand, it seems that the proper place of the RTC in the framework of the general theory of the integrable 
systems has not been adequately understood yet. Since the seminal papers of Kyoto school [0 it is evident that any integrable 
system can be treated as a particular reduction of the (multi-component) Kadomtsev-Petviashvili hierarchy. This statement 
is certainly true for the RTC too. Using the powerful machinery of fl|, it is possible to describe wide classes of solutions to 
the RTC and reveal many interesting algebraic properties of the system. Another important lesson of the same philosophy 
teaches us that many would-be different systems have the common structures behind if considered on the level of fundamental 
notions like Baker-Akhiezer functions, Grassmannian manifolds etc. In this paper we want to advocate these points more (in 
somehow reversed order). Indeed, we shall show that the RTC has the deep connection with many well known integrable 
systems. Namely, the special case of the AKNS hierarchy Q as well as the discrete non-linear Schrodingcr equation (of type 2) 
Q ^ are nothing but different faces of the RTC. Some "novel hierarchy" considered recently |jj is equivalent to RTC also; we 
shall see that all these systems belong to the same class of equivalence: the corresponding L-operators are connected with each 
other by the suitable gauge transformations. Another interesting observation is that the forced (i.e. semi-infinite) variant of 
the RTC is described by the unitary matrix model which, in turn, naturally arises in the context of the two-dimensional gravity 
|lo[ pj| . Actually, the later case was our starting point. Namely, we have started from the integrable system (with the infinite 
set of the integrable flows) J12] ] inspired by the unitary matrix model. The equivalence with the discrete AKNS hierarchy was 
established there. Here we show that this system is nothing but the generalized version of the RTC hierarchy. 

The connection with matrix models ^] turns out to be important by some reasons. First of all, it allows one to apply effectively 
the orthogonal polynomials technique, which considerably simplifies all the derivations and proofs. Indeed, in the systems of 
such a type there are two different (but essentially equivalent) Lax representations - those given by 2 X 2 matrices depending on 
the spectral parameter and given by infinite matrices (in spirit of |l3|). We mostly base on the second representation, which is 
effectively treated by orthogonal polynomials methods (since they solve the corresponding Riemann-Hilbert problem - see [pl|). 
As a by-product of our consideration, we obtain some new results specifically interesting for the orthogonal polynomials theory. 
Second, it was shown in Jl2| that unitary matrix model is equivalent to the two-component Volterra system, thus establishing 
the connection between this latter and the RTC. On the other hand, Volterra systems are well investigated and, therefore, this 
last connection might help to get more about the RTC system from this different point of view. 

There is another even more useful representation of the RTC. Namely, we demonstrate in this paper that the RTC is a 
reduction of the two-dimensional (subholonomic^]) Toda lattice hierarchy (2DTL). On this level, the equivalence of the above 
mentioned systems is especially transparent. Looking at the RTC as a reduction of the 2DTL also allows us to give some r- 
function interpretation of the RTC and realize it in terms of fermionic correlators. This finally implies the bilinear identities like 
those proposed in Q and the theory group interpretation. Let us also stress that the r-function (Grassmannian) interpretation 
being absolutely algebraic can be easily deformed to the quantum case (along the line of the papers JwJ,J^]) Q Furthermore, 
using group-theoretical approach one can deal with different RTC's - that is, the general RTC, forced RTC and RTC molecule 
hierarchy - from the algebraic point of view. Their algebraic properties are very interesting. In particular, the latter hierarchy is 
related to a class of co-representations of U q (SL(N)) in the quasi-classical limit. Put it differently, in the Hamiltonian approach 

Let us stress that the correspondence of the RTC and a unitary matrix model is not so surprising: it is well-known that 
the Hermitean matrix model describes the TC hierarchy . On the other hand, the RTC is usually considered as a group 
generalization of the TC whatever it means. Therefore, it is natural to identify the RTC with some unitary matrix model. 
2 See 

3 Quantization of the RTC was considered in W], where the Sklyanin's scheme of the separation of variables was applied. 
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the relativistic Toda molecule is given as the integrable system on some Poisson leaves on SL(N) (see Jit| for mathematical 
background). This is why we consider this approach as a very promising one. 

This paper is organized as follows. In section 2 we discuss the simple spectral problem (three-term recurrent relation) for 
the RTC which is natural generalization of the corre sponding problem for the usual Toda Chain (TC). Then, we introduce two 
different flows which lead to the same equation (1.1 ) and interpret them in the context of the results by Suris |^]. These results, 
we guess, arc not new and can be partly extracted from preceeding papers (see, for example, ji| ^|). The main aim of section 2 
is, however, to establish the key points which allow us to connect all the integrable structures mentioned above. 

Section 3 contains the derivation of the both RTC spectral problem (proposed in Section 2) and the evolution equations 
from the (bi)orthogonality conditions implied by the unitary matrix model. As a by-product of the orthogonality conditions, 
remarkably, we automatically get the spectral problem in such a form that describes the connection of the Baker- Akhiczcr 
functions of the two-dimensional Toda Lattice hierarchy. We show that the solution to the RTC hierarchy satisfies simultaneously 
the first 2DTL equation with respect two the flows introduced in the previous section. The connection with the two-component 
Volterra hierarchy as well as with the AKNS system is also discussed. As a consequence, we get the direct correspondence 
between the RTC, a particular discrete version of the non-linear Schrodinger equation, and the "novel" integrable system |J. 
At the end of the section, we reveal the existence of two Lax operators for the RTC which naturally arise from an immediate 
generalization of the recurrent relations inspired by the unitary matrix model. 

In section 4 we outline the essential ingredients of the 2DTL theory [ [l3| and explicitly describe the whole RTC hierarchy 
as a particular reduction of 2DTL. 

The forced RTC hierarchy and its finite analog, the relativistic Toda molecule, are discussed in sections 5 and 6 respectively. 
We give the group interpretation of the later hierarchy. We also discuss the fermionic representations for the r-functions of the 
corresponding hierarchies and manifestly indicate the reductions from the 2DTL which describe the hierarchies in terms of the 
Grassmannian. We present also the explicit solution to the forced RTC in the determinant form. The fermionic representation 
received allows one to make some connection with the results of jij , but we postpone the discussion of this point till a separate 
publication (see, however, section 9). 

In section 7 we describe the simple approach to discrete evolutions of the RTC which is based on the notion of the Darboux- 
Backlund transformations and can be considered as a natural generalization of the corresponding notion in the usual Toda chain 
theory. The continuum limit as well as the limit to the Toda chain is also discussed here. We consider also some degeneration 
of the general Darboux-Backlund transformation leading to simplified discrete equations which can be treated as a kind of 
"modified Toda chain" (these equations have been proposed for the first time in |J but we should stress that they have nothing 
to do with the usual discrete-time Toda chain since the Lax operators of these two systems are different). 

Section 8 is devoted to some applications to the theory of biorthogonal polynomials. Indeed, we discuss the (relativistic) 
orthogonal polynomials which satisfy the three-term recurrent relation specific for the RTC and look at some very specific 
solutions to these relations. One of these solutions describes finite sets of the polynomials corresponding to the Toda molecule, 
the other one is a special limit of the Askey polynomials ]l8| . In the section, we discuss the orthogonality measures and some 
additional orthogonality relations which can be manifestly treated in the indicated simple cases. 

In concluding remarks we briefly discuss more general (and, in a sense, more natural) description of the RTC hierarchy 
which will be considered in detail in the forthcoming publication Jis| ] . 

In Appendices we give the proof of the some important assertions used in the main body of the paper. In particular, in 
Appendix A there were obtained the evolution equations starting from the orthogonality conditions determining the unitary 
matrix model, while in Appendix B the spectral problem naturally arising in the framework of orthogonal polynomials is 
transformed to that of the 2DTL hierarchy. 



2 Lax representation for RTC 

In this section we describe the Lax representation for the standard RTC equation. The usual procedure to obtain integrable 
non-linear equations consists of the two essential steps: 

i) To find appropriate spectral problem for the Bakcr-Akhiczcr function(s). 

ii) To define the proper evolution of this function with respect to isospectral deformations. 
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2.1 Lax representation by three-term recurrent relation 

In the theory of the usual Toda chain the first step implies the discretized version of the Schrodinger equation (see ||, for 
example). In order to get the relativistic extension of the Toda equations, one should consider the following "unusual" spectral 
problem 

*n+l(«) + a„0„(z) = «{* n («) + &n0 n _i(z)} - n e Z (2.1) 
representing a particular discrete Lax operator acting on the Baker- Akhiezer function <J> n (z). This is a simple three-term 
recurrent relation (similar to those for the Toda and Volterra chains) but with "unusual" spectral dependence. 



As for the second step, one should note that there exist two distinct integrable flows leading to the same equation (1.1) 



As we shall see below, the spectral problem (2.1) can be naturally incorporated into the theory of two-dimensional Toda lattice 
(2DTL) which describes the evolution with respect to two (infinite) sets of times (ii, •••) i {t— lit— 2i •••) (positive and 
negative times, in accordance with Jli|). In section 4, we shall derive the general evolution equations satisfied by <E> n . Here we 
describe the two particular flows (at the mome nt, w e deal with them "by hands", i.e. introducing the corresponding Lax pairs 
by a guess) which lead to the RTC equations (1.1). It turns out (as usual) that the simplest equations are to be associated 
with the first times of this hierarchy ti and t—i. The most simple evolution equation is that with respect to the first negative 
time and has the form 

Rn^n-l (2.2) 



dt-i 

with some (yet unknown) R n . Remarkably enough that the form of this equation exactly coincides with the corresponding 
equation for the non-relativistic Toda chain. 

The compatibility condition determines Rn in terms of a n and b n : 

b n 
a„ 



Rn 



(2.3) 



and leads to the following equations of motion: 



da n 
<9t_l 
db„ 



dt- 



«n+l 

O.n + 1 

_L_ 

O-n — 1 &n 



i 



In order to get (1.1), we should identify 



a n = exp(-ep n ) 
-<e 2 exp(g n - q n ~i - ep n ) 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



(2.8) 



Note that in this parameterization the "Hamiltonian" R n in (2.3) depends only on coordinates q n 's 

Rn = exp(g n - <? n -i) 

Such a form is typical for integrable systems (for example, the first flow of the TC has the same structure). 

In order to preserve the proper limit to the TC equations as t — » 0, one should perform the rescaling of time in (2.4), (2.5) 

*_1 ^ K<0*-1 ( 2 -9) 
where the function v(e) is to behave as 1/e in the limit of small e. Then, from ( |2,4| ), ( |2.5[ )p| 

dg n 



dt- 

dpn 

dt- 



After choosing 



-i^(e)|l + e 2 exp(q n+ i - g n )| cxp(ep n ) + v(e)(l + e 2 ) 
/(e)|exp(<?„ - q n -i + ep n -i) - exp(g n+ i - g n + ep n ) j 



(2.10) 
(2.11) 
(2.12) 



4 Strictly speaking, ( |2.5| ) gives the equation g^ n 



9q n -l 



: Qn — Qn—l, where Q n is the exponential term in (2.10). Thus 



one can only write 



-u(e){l + e 2 exp(<j n+ i - q n )} exp(ep„) + constant 



dt-! 

For the fast-decreasing solutions, constant = u(e)(l + e 2 ). This choice of the constant nicely agrees with the non-relativistic 
limit e — > 0: in this limit g^ n = p n - 



G 



these equations trivially lead to the RTC equation 

d 2 q. 



d In (-. . dq n+ i\ ( dq n \ 



dq„ \ exp(<j n+ i - q n ) 



1 + e 



at- 



l + e 



9t_l / 1 + e 2 exp(q n+1 - <j n ) 
9?n-i \ exp(g n - q n -i) 



(2.13) 



dt-i ) 1 + e 2 exp(g n - q n -i) 



Another form of the same equation is given by the replace 

?n = — = exp(ep n ) 



(2.14) 



fl?) fl?) — 1 



Then, from (p!4), O ) | 



-e 2 exp(q n - q„-i + ep n -l) 



■ £n(»7n+l - »7n) 



(2.15) 



: Vn(Vn+l - Vn-1 + fw-1 - Cm) 



As we noted already, the evolution (2.2), which leads to the RTC equations is not the unique one. The other possible choice 
leading to the same equations is 

^ =-&n($n-**n-l) (2.16) 
oti 

The compatibility condition of ( j2 . l| ) and (2.16) gives the equations 



-Qn(fen + 1 — b n ) 



TT— = -»n(«n+l ~ O-ra-l + a n-l — "n) 
Ot l 



(2.17) 



(2.18) 



which are exactly of the form (2.15). In terms of (q n , p n ) we have the following equations (with the same rescaling of time as 
in (E3)): 

,),.. f i 

(2.19) 



= -K e )|l + ex P(<?" - qn-l)\ cxp(-ep„) + i/(e)(l + e 2 ) 
"7~~~~ = e^(e)|exp(<j, l+ i - q„ - ep n+1 ) - exp(q n - g„_i - ep„)| 



(2.20) 



slightly different from (2.1C), (2.11). Nevertheless, the second order equation for q n is exactly (2.13). 



2.2 2x2 matrix Lax representation 

The same RTC equation can be obtained from the matrix Lax operator depending on the spectral parameter (generalizing 
the Lax operator for the TC Q). Then the RTC arises as the compatibility condition for the following 2x2 matrix equations: 

inVn = V'n+l , —^7- = A n 1p n (2-21) 

at 



where 

r(S) 



£exp(ep„)-( 1 eexp(cj„) 
-ecxp(-q n + ep n ) 
/ <e 2 exp(g n - q n -i + epn-l) 




Using (2.22), one can re-write the matrix spectral problem as the recurrent relation for i/>„ 



(2.22) 
(2.23) 



< 1 cxp(-q„_i + epn-i) 
ictral problem as the recur 
V>L+i(C) = {C exp(ep„) - C^nHt) - e 2 exp(<z„ - fo-i + ep„-l)^ 1 2i(0 (2.24) 

Then, putting 

^'(OhKMO ( 2 - 25 ) 

5 In the paper [gj d n = £ n , c„ = — »? n +l >t = — i— 1 and e = 1. 
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where 



one sets 



n+l 



: £exp(-ep n ) 



*n+l(C) + exp(-ep„)*„(C) = C 2 {*n(C) - € 2 exp(g n - g„_i - ep n )$„_i(C)} 
This is exactly recurrent relation (2.1) with a n , b n given by (2.6), ( |2 . t[ ) after the identification 

z = e 

In the same manner, the evolution equation in terms of <!>„ takes the form 

—— = e exp(q n - <j n _i)$„_i 
at 



Thus, the evolution determined by (2.23) is associated with —t—\ -flow in our approach (see (2.2), (2.8)). 
Obviously, we are able to construct another A-operator which generates a new integrable flow: 

/ -e(cxp(q n - ep n ) \ 



An 



<exp(-g n _i) C 2 + e 2 cxp(q n - q n ^ 1 - ep n ) - 1 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



(2.30) 



It is easy to see that (2.30) gives the evolution equivalent to (2.16) on the level of the Baker-Akhiezer function <E>„. 



To conclude this section, we should remark that L-operator (2.22), which determines the RTC is not unique; moreover, 
it is not the simplest one. Indeed, in the next section we shall see that there exists the whole family of the gauge equivalent 
operators, which contains more "natural" ones and includes, in particular, the well known operator generating the AKNS 
hierarchy. From general point of view, the whole RTC hierarchy is nothing but AKNS and vice versa. 



3 RTC and unitary matrix model, AKNS, etc. 

Now we are going to describe the generalized RTC hierarchy as well as its connection with some other integrable systems. We 
start our investigation from the framework of orthogonal polynomials. The advantage of it is that one does not need to guess 
Lax pair, but instead can get it automatically using the formalism developed in [0, |io| . From the point of view of the 
integrable systems, this method is nothing but the Riemann-Hilbert problem Jll], |20|. The orthogonal polynomials give the 
solutions only for hierarchies of the forced type (see the proper definition below). Fortunately, the hierarchies we consider are 
"local" ones (one-dimensional subfamilies of the subholonomic systems [^|) either can be reduced to these, and, therefore, Lax 
representations obtained from the orthogonal polynomials formalism can be continued to the whole hierarchies pp| . From now 
on, we mostly arbitrarily turn from the forced hierarchies to the general ones and back. However, right in the RTC case, this 
point still deserves some comments which we postpone till section 5. 



3.1 Unitary matrix model 

It is well-known that the partition function r n of the unitary one-matrix model can be presented as a product of norms of the 
biorthogonal polynomial system |2l[ . Namely, let us introduce a scalar product of the form ^| 



< A,B >= I ^^exp( V(t m z m -t- m z- m )\A(z)B(-) 



(3.1) 

m>0 

where the integration measure is not fixed but to be chosen so that the integral in (3.1) is well defined; for example, if dfi(z) 
can be presented in the form dz exp(^ c k z k ) then everything reduces to the usual measure dz with the proper shifts of times 
ifc , t_ fc. Let us define the system of polynomials biorthogonal with respect to this scalar product 

< * n , >= h n S nk (3.2) 
Then, the partition function t„ of the unitary matrix model is equal to the product of h n 's: 



n-l 

T„ = 

k=0 



Yl h k , r = 1 (3.3) 



6 The signs of positive and negative times are defined in this way to get the exact correspondence with the times introduced 
in El. 
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the integration contour in ( fc.l| ) being the unit circle^] (and, therefore, h„'s are real for the real times t^, t_k's). The polynomials 
are normalized as follows (we should stress that superscript '*' does not mean the complex conjugation): 



*n(«) 



■ + S„ 



*n(*) 



z n + ... + s* 



(3.4) 



Now it is easy to show that these polynomials satisfy the following recurrent relations, which gives the discrete Lax operator 

*n+x(«) = z$ n (z) + Snz"®^- 1 ) 



(3.5) 



and 



(3.6) 



These recurrent relations can be written in several equivalent forms. First, it can be presented in the form analogous to 
( |2.l| ), i.e. equivalent to the spectral problem for the RTC: 

*n+l - ^-(1 - S„_lS£_i)*„_l) (3.7) 

->n-l I On-1 > 



4>*. 



s ... 



<I>" 



Mi-£.-i3U)*£-i 



From the first relation and using (2.27), one can immediately read off 

^ n 



= -cxp(-ep n ) 
<? exp(<j n - q n -i) 



(3.8) 



(3.9) 



h„ = 7(-f 2 ) n exp(g n ) 



(3.10) 

where 7 is a constant which does not depend on n. Thus, the orthogonality conditions (3.2) lead exactly to the spectral problem 
for the RTC. We should stress that equations (3.7), (3.8) can be derived from system (3.5). Usually, the three-term recurrent 
relations can be transformed to the system of two equations in many different ways (see below). The main feature of (3.5), 
which distinguishes it from the other possible choices is the remarkable fact: this system describes the connection between the 
pair of the Baker-Akhiezer functions arising in the context of 2DTL. We shall discuss this point in section 4. 

Some remarks are in order now. In all formulas above the discrete index n runs over the non-negative integers. One 
can trivially extend all relations (3.5 )-( 3. 10 ) (and all the relations below) to arbitrary n 6 TL . Then, the polynomial case 
corresponds to the conditions 

- r ^— = , n < (3.11) 

(3.12) 

(3.13) 



or, (see (3.9)) 



1-1 



with (formal) ordering g_i < g_2 < • • ■ • Equivalently, the analytical continuation of (j?^) written in the form 

h Tn+1 

tin = 



reduces to the polynomial case by imposing the conditions 

r„H0, n < (3.14) 

This is exactly what we call the forced integrable systems pfj| . Our way of doing is to derive all the formulas in the forced (= 
polynomial) case and then to extend them to arbitrary values of the discrete indices. 

Using the orthogonal conditions, it is also possible to obtain the equations which describe the time dependence of <J? n , . 
Namely, differentiating (3.2) with respect to times t\ , t—i gives the following evolution equations (for derivation see Appendix 
A): 

d$„ S n h„ ,^ ^ 

(3.15) 



Sn~l hn — 1 



dt-1 
^ n 



h n -l 



($ n - z* n _i) 

-*7!-l 
<f"*_1 



(3.16) 
(3.17) 



7 We do not restrict ourselves to this integration contour throughout this paper. 
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9$* 



S* h n 



71—1/ 



(3.18) 



9t_l 'S'n-l ^i-l 

(see general evolution equations with respect to higher flows in the next section). The compatibility conditions give the following 
nonlinear evolution equations P: 



As a consequence, in the polynomial case, 



9S n h n +l 

- b n +l 
oil h n 


(3.19) 


riQ h , ^ 
OJ n c ft-n+1 

Ot-1 /in 


(3.20) 


95* _ + /ln+1 
9*1 U_1 /in 


(3.21) 


95* c * /i„+i 

Ot-l /in 


(3.22) 


9/ln * 


(3.23) 


9/l n i 

„ — S n —±S n h n 


(3.24) 



These are exactly relativistic Toda equations written in somewhat different form. Indeed, from (3.23), (3.19) and (3.21) one 
gets 

9 2 



Using (3.23) again and ( fc.rj ) 
Substitution of SnS^_ 2 to (3.25) gives ^ 



^ log /In 



hn + 1 „ „* /lr, 

— h OnJ„_2T 

/In "ti- 



/ 5 

V9iT log ' ! " 



(- 

V9*i 



log/i n _i = S„5. 



V hn-lJ 



(3.25) 



(3.26) 



/^n^ 



92 , I, 

^loghn 



(a7T log/l " 



f- 

V9ti 



log ^ + 1 



1 - 



(3.27) 



/in 



On the other hand, the RTC is a particular case of the 2DTL hierarchy. Indeed, let us introduce the key objects in the theory 
of integrable systems - the r-function as it is defined in (3.13). Then, with the help of (3.19)-(3.2^ ), one can show that the 
functions 



satisfy the first equation of the 2DTL: 

dt 1 d t _ 1 logTn 



(3.28) 



Therefore, it is natural to assume that the higher flows generate the whole set of non-linear equations of the 2DTL in spirit of 
|l3[ . In the next section, we shall see that this is indeed the case. On the other hand, it turns out that the integrable system 
dete rm ine d by (3.7), (3.8) is highly degenerate as compared with the general 2DTL. Indeed, the solutions to the equations 
( |3.7| ), ( [3.8! ) (treated for a moment as two separate equations) correspond to the four Baker-Akhiezer functions of the 2DTL. In 
the general theory [jl^, these functions are linearly independent. In the RTC case, they are dependent due to (3.5); this is the 
origin of the degeneracy. Actually, this degeneracy is responsible for th e appearance of additional non-linear equations, which 
are absent in typical situation. For example, these are equations (3.27), which contain the derivatives of log/i n with respect 
to ti (or *_i) only. These additional equations being consistent with the whole 2DTL hierarchy specialize the reduction of 
the 2DTL. Therefore, we should treat the RTC hierarchy as a special reduction of an "abstract" 2DTL (see section 4). |^| 

Let us also note that, using (3.10), one can get from equation (3.27) the equation distinguished from (1.1) by the linear (in 
time) shift of the variable g„ — > q„ — —ti. This means that one should careful when express q„ through T-functions. 

This completes the derivation of the RTC from the unitary matrix model. Now we consider the connection of the RTC 
with the integrable systems mentioned in the Introduction. 

8 In the po lynomial case, one can get these equations directly from (3.1!: )-(3.18 ) considering terms of the order ~ z°; see 
definition (Q. 

9 The same equation holds for t_i-flow. 

10 Le t us s tress th at th e term "reduction" does not obligatory mean the polynomial case since non-polynomial solutions to 
( |3, ?] ), ( |3,8| ), ( |3. 15| )- (3.18) are also reduced in the above mentioned sense. 
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3.2 Two-component Volterra 



Recursion relations (3,7)-( p.8| ) can be presented in the form establishing their connection with the two-component Volterra 
hierarchy |l2| (it should be understood that, generally, <E> n and <E>* in these relations are not polynomials). For doing this, we 
introduce the two-component functions 



f i n ±) {z) = z-^+ 1 ^ 1 (z)±z^~ 1 ^ n _ 1 {z- 1 ) 
These functions satisfy the following orthogonality conditions: 

< /i ±) ,4 ±) >= [2T(S n -l+Stl)] hn-lSnfc 

< fk +) ,ft } >= (Sn-i - S^_ 1 )/ l „_i5 nfc 



(3.29) 



(3.30) 



These functions also satisfy the recurrent relations which can be written in the compact matrix form after introducing the 
two-component column 



f(+) 

J n 
.1 n 



Then, one gets the recurrent relations giving the two-component Volterra Lax operator 

(z3 + z~i)f n = + V„fn-1 

(Z2 - Z~2)f n = a 1 f n + l + Vnfn-l 

where 



cr;'s are Pauli matrices and 



Wn = h n -i \ ; det Wn = 2hnh n -i 

V -(1 + Sn-l) 1 + S£_1 / 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



Three forms of the recurrent relations (3.5), (3.7) and (3.32) explicitly establish the connection between the unitary matrix 
model, RTC and two-component Volterra system (certainly, there still remains to compare the evolutions, see the next section). 



3.3 RTC versus AKNS and "novel" hierarchies 

Now let us demonstrate the correspondence between RTC and AKNS system. We have already seen that the orthogonality 
conditions naturally lead to the 2x2 formulation of the problem generated by the unitary matrix model: 



L (u) 



a 



K+i 



L<°> 



z n S n 
z- 1 



<1>„ 
<1>* 



S n S n _ 1 



d 



-i 



-n 



Z Sn 



Sn—lS n 



(Equations (3.3C), (3.37) follow from ( 3. 15 )-( 3. 18 ) and the original spectral problem ( p\5| )). Put 

<T>* = z -"/2+l/4 F * 



Then the spectral problem (3.35) can be rewritten in the matrix form 

T (AKNS) -77 -77 .77 



F* 



where 



r (AKNS) 



s* C 1 



,1/2 



(3.35) 

(3.36) 
(3.37) 

(3.38) 

(3.39) 
(3.40) 
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This is the Lax operator for the discrete AKNS j7J. Obviously the evolution equations ( |3.36[ ) , (3.37) can be written in terms of 
F„, F* as 



9*1 ~ " >n 



at. 



( ~C 2 C^n-l 



Note that after introducing the trivial flow 



,(0) 



n o -i ; 

we get the difference non-linear Schrodinger system (DNLS) Ti| (see also H) generated by the "mixed" flow 



d:F " _ 4(0)77 
OT 



(3.41) 
(3.42) 

(3.43) 



dr -\dto at-! atJ^ n - (An +An An Wn 
( \ + S n S* n _ x -C % C _1 5„-i-C5„ \ 

V C 1 S n — C^n— 1 — ^ — S n —lS^ + C 2 / 



(3.44) 



Indeed, from the compatibility conditions for (3.401, (3.44) or, equivalently, just from ( J3. 19| )-( 3 . 2C ) (along with the trivial 
evolution dt S n = 2S n , dt S n = —25* ) one gets the discrete version of the nonlinear Schrodinger equation: 



BS. 



— ~(Sn+l — 2S n + S n -l) + SnS n (Sn+1 + Sn-l) 



(3.45) 



Note also that the "novel" hierarchy of |9j is equivalent to the RTC (and, therefore, to the AKNS hierarchy) as well. Namely, 
the Lax operator in |1, i.e. 



Z -\- Un'Vn t^n 



b (1) 

"n+1 



1 



6 (2) 
"n+1 



defines the recurrent relation of the form (3.7): 



^ll - f u n v n + *L.) tf> = Z ( & - -^^O 

\ Un-XJ \ Un-1 / 



thus revealing the connection with the RTC. Comparing (3.7) and (3.47) leads to the identification 



Un — Snhn , 



where h n 's satisfy (^). Moreover, from ( fc.46| ) and (Ej) it is easy to sec that 



and, therefore, L n can be obtained from L„ AKN by the discrete gauge transformation: 



T —TT r(AKNS) rr -l 
L<n — Un + l^n u n 



(for z = f 2 ) where 



U n = z n/a-i/4 



1 





S*n-1 






h n 



- U n +1 — U n V„ 



du„ 



Evolution equations (3.19)- (3.24) in terms of new variables (3.45) have the form 



U n -1 



-V n -1 + UnV n 



dt-l 1 + Un-lVn 



dh " dt-l 1+UnVn + l 

and easily reproduce the usual AKNS equations in the continuum limit since 

(9t —dt-i- 9t-l)«n = - - 2u n + U n -l) + (" n _l + «„) V n + ■ 

(9t - 8t 1 ~ dt_ 1 S )v n = (vn + 1 - 2V„ + Un-l) - (l>„ + + + ■ ■ 



(3.46) 



(3.47) 



(3.48) 



(3.49) 



(3.50) 



(3.51) 



(3.52) 



(3.53) 
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We conclude with the remark that the operator L„ in (2.22) is also gauge equivalent to L$j AKN 



where 



r(S) _ Tj , . r(AKNS) fr l 



( tl)l 

S n —i 

2n-l ,-1/2 



\ Sn-lhr, 



\ 



(3.54) 



(3.55) 



3.4 Non-local Lax representation 

At the end of the section let us note that there is another form of the recurrent relations which is non-local (i.e. contains all 
the functions with smaller indices) but instead expresses <J» n (z) through themselves. This form of the spectral problem will be 
crucial for dealing with the RTC as a particular reduction of the 2DTL. Let us introduce the normalized functions 



V n (z) = <S> n (z) 

such that, for example, in the polynomial case 



(3.56) 



(V n ,V£)=5 nk (3.57) 
Then, from ( |3.7| ), ( |3.8| ) one can show that in the forced and fast-decreasing cases some proper solutions (see the discussion in 
Appendix B) satisfy the equations 



zr n (z) = r n+1 (z)-Snh n V -±±v h ( z ) = c nk v k (z) 



- 1 P*(^ 1 ) = ^±ip* +1 (z- 1 )-5: Sk-iVti*- 1 ) = CknVU*- 1 ) 



(3.58) 



/in 



Note that this expression is correct for general (non-polynomial) Vn and T 3 * provided the sums run over all integer k. In the 
polynomial case, the sums automatically run over only non-negative k. Let us also stress the natural appearance of variables 
(3.46) in the Lax operator C nk . The last representation of the spectral problem will be useful in the next section to determine 



the general evolution of the system. Indeed, these relations manifestly describe the embedding of the RTC into the 2DTL 
(l3l [13], which is given essentially by two Lax operators (C and C). 



4 RTC as reduction of 2DTL 
4.1 Hierarchy evolution 

In order to determine the whole set of the evolution equations, one can use different tricks. The simplest one is to use the 
orthogonal polynomial technique and then to continue the result for the "general" hierarchy. We use this way at the end of the 
section. Now let us note that there is also some less direct way to obtain the evolution suitable even in the non-polynomial 



case without continuation. Namely, one can use embedding (3.58) of the system into the 2DTL, making use of the standard 
evolution of this latter |l3| . Let us briefly describe the formalism of the 2DTL following . In their framework, one introduces 
two different Baker-Akhiezcr (BA) ZxZ matrices W and W. (In the original paper JlS] , these matrices were denoted as W^°°^ 
and ; the superscripts correspond to the marked points on the Riemann surface which specify the solution to the 2DTL, 

i.e. the point of the infinite dimensional Grassmannian) . These matrices satisfy the linear system: 

i) the matrix version of the spectral problem: 

OV = WA , CW = WA- 1 (4.1) 

ii) the matrix version of the evolution equations: 

_ (4-2) 
^- = (£ m )_W, ^L = (C m )-W, m = l,2,... 
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where Z X Z matrices £ and £ have (by definition) the following structure: 

£ = ^diag[6 l (s)]A l ; 6i(s) = l 

i<l 

(4.3) 

C = ^ diag{ Cl (s)]A l ■ c_i(s) £ 
i>-i 

Here diag[fc;(s)] denotes an infinite diagonal matrix diag(. . . 6j(— 1), 6^ (0) , fci(l), . . .) ; A is the shift matrix with the 
elements A nk = S nk _i and for arbitrary infinite matrix A = ^? diag[di(s)]A l we set 

s ^diagMs^A 1 , (A)_ = ^dia«[oi(a)]A' (4 4) 

i.e. is the upper triangular part of the matrix A (including the main diagonal) while (A)_ is strictly the lower triangular 

part. 



Note that (4.3) can be written in components as 

C-nk = <Wi,fc + b k- n (n)0(n - k) , C„ k = c_i (n)<5 n _i ifc + c k _ n {n)0{k - n) ; n,keZ (4.5) 
The compatibility conditions imposed on ( |4.l| ),( |4~^ ) give rise to the infinite set (hierarchy) of nonlinear equations for the 
operators C, C or, equivalently, for the coefficients b m (n), c m (n). This is what is called 2DTL hierarchy. 

On the level of nonlinear equations, one does not need the information of the structure of BA matrices. However, in order 
to get the touch with the polynomials, this information is essential. It was proved in Jl3| that the linear system ( |4,l| ), ( |4.2[ ) is 
resolved by the following BA matrices: 

W = Wexpi jTt m A m 1 , W = Wexpi ^ t_ m A~ m 1 (4.6) 

^ m— 1 ) \ m=l ) 

where W, W can be presented in the form: 

oo oo 

W = V" diag^^A^ , W" 1 = Y A_i diag[<(s + 1)] , w (s) = w%(s) = 1 (4.7) 

i=0 i=0 

oo oo 

W = V" diaglTuiCs^A* , W - 1 = V" A'diagpo^s + 1)] (4.8) 

i=0 i=0 
Now let us introduce the BA functions as follows: 

oo 

w n (z) = V" W nk z k = z n ^^ Wi (n)z- 1 (4.9) 
fcsZ 1=0 

OC 

w n {z) = ^2 W„ fc z fc = z n '^w i (n)z i ( 410 ) 

fceZ l =» 
At the same time, we define the adjoint BA functions as 

oo 

ZW* + 1 (z) = Y W kn Z ~ k S w n( z )= z ~ n ^2 w i( n ) z ~ i (4.11) 

fceZ 

OO 

zw* +1 (z) = ^ W fcn 1 z~ fc ; w*(z) = z- n ^^w*{n)z l ( 412 ) 

feeZ '=° 
Using ( |4.l[ ), (^3), it is easy to show that the BA functions satisfy the linear equations: 

C„ k w k (z) = zw n (z) , C nk w k (z) = z^ 1 w n (z) (4.13) 



9w n (z) dw n (z) 

—^7 = -[(£ )-]nk«'k(*) , — = [(£ ) + ] nfc u> fe (z) 

~ l(£ )-Jn*">fcO*) i ^7 = -[(£ ) + J nfc u> fc (2:) , m = l,2, 



(4.14) 



Corresponding system for the adjoint BA functions is: 

C kn wl +1 (z) = zw^ +1 (z) , C kn w k * +1 (z) = z _1 ijr* +1 («) (4.15) 



Sw*,, ( z ) dw*,-,(z) 

" + = [(O-Un^ + lW , " + U = -[(£ m ) + ] fc „^* +1 ( 2 ) 

Otjn Otm 

9 W n-Ll( z ) —m 9w* +1 _ m _ 

— ^ = -[(-£ )-]nfcWfe+i(z) , — = [(£ )+]knW k+1 (z) , m = l,2, 

Ot— m OZ — 



(4.16) 
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4.2 Reduction from 2DTL 



Now let us return to the case of RTC "in general situation" , i.e. with the conditions S n = 1, S£ 
polynomial case being removed away. From (3.56), one gets two matrices 



S n +l,k - T^SnSk-xOin- k) , fe,neZ 



1 n < restricting onto the 



(4.17) 



-<5n-l,/s ~ S n -lS£9(k - n) 



k,n 6 Z 



(4.18) 



which have exactly the form ([0 ) . One can consider the spectral problem (4.13) for these particular operators. It is easy to 
prove that every solution to (4.1J), (4.17), ( |4,18| ) satisfies the recurrent equation (3.7). On the other hand, ( |j.7j ) has two linear 
independent solutions with asymptotics (we are using the identification (3.5C)) 

V£\z) = h n z n (l + 0(z)^ z^O 



(4.19) 



(4.20) 



which are precisely the same as asymptotics of (4.9) and (4.1C) respectively. Thus, the solutions to (3.7) should be identified 
with the corresponding BA functions w n (z) and w n (z): 



w n {z) = V^\z) 
w n {z) = Vi 2) (z) 



(4.21) 
(4.22) 

Similarly, one can consider the spectral problem ( |4.1E| ), (1.17), (4.18 ). T hen, t he c orresponding w^(z) and w*(z) satisfy the 

(4.23) 



same three-term recurrent equation as V^iz 1 ) = <!>*(£ 1 )/h n (see (3.56) and (3.8)), i.e. 



h n -\ 



' n+l 



s 



12_-p* 



This equation obviously has two independent solutions with asymptotics 



p: (1) ( 2 - 1 ) = 2 -"(l + o(i)) z^oo 
V* (2) (z- 1 ) = l- z - n (l + 0{z)^ z^O 



(4.24) 
(4.25) 



which are exactly of the form (4.11) and (1.12). Thus, the adjoint BA functions should be identified with solutions to (4.23) as 
follows: 

w*(z) = z- 1 V^\{z~ 1 ) (4.26) 
w:(z) = z- 1 V* i ? 1 (z- 1 ) (4.27) 



To conclude, we see that equation (3.5£ ) (see also (4.26) below) for "P n (z), V^(z 1 ) is exactly the same as that for w n +i(z), 
w* +1 (z) (compare with (4.14), (4.16)). 



4.3 Evolution and orthogonal polynomials 



Let us now demonstrate how one can obtain evolution of the RTC using the technique of the orthogonal p olyno mials. Indeed, 
differentiating the orthogonality conditions with respect to arbitrary times, one can obtain with the help of (3.58) the evolution 
of polynomials V n and (for details see Appendix A) : 



&Pn 
dt m 

dT n 

dt-m 

^ ' n 
dt m 

8V* 

dt-m 

, — )nn ll "n 
Ot m 



[(C m )-] nk V k 

[(c m )+] kn vi 

dh„ 



(4.28) 



9t_ 
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where by definition for given matrix C n ki 

[(C) + } nk = C nk d(k-n) , (C-) nk = C nk e(n-k-l) 
i.e. (C)-|- is the upper triangular part of C while (C)_ is strictly the lower triangular part of C 



(4.29) 



Let us note that this evolution (of the BA functio ns) h as the same form as (4.14), (4.16). However, the latter one was 
determined on the infinite matrices, while the evolution ( 4.28 ) - on the quarter part of these matrices (this is the forced hierarchy, 
we will return to this point later). In order to get the hierarchy of non-linear equations, one should use the compatibility 
conditions of the linear systems. The hierarchies obtained in this way from the two considered evolutions do not coincide 
exactly although differing only by inessential constants. Now we consider in more details the evolution with respect to the first 
times to observe this phenomenon. 



Using the general formulas (4.14), (4.16) one gets the simplest evolution equations 

71 — 1 ~. 



dti 



<9$n 
8t-l 

9K = 



Sn- 
h n 
hn-1 



-(*„ 



h k 



-<£V 



-l h n -i 

<!> n _l = (1 - Sn-lS*_l)*n-l 



*n-l = -(l- 



dt-_ 



~ S n hn 



n-1 

E 



Sn—lS^ 



-l)*n- 



hk 



(*;, 



n — 1 / 



The compatibility conditions give the following nonlinear equations: 



8Sn _ 

at! 


c hn + 1 
h n 


f aS n 


dSn 

at-! 


O hn+1 
On-1 , 

hn 


+ PS n 


as* _ 
att 


°n-l , 

h n 


- a S n 


9S+ 
dt-x 


a* h n +i 

lln 


-I3S* 


dh n 
~dh~ ~ 


— {S n S n _ 1 H 


- l)h n 


dh n 

dt-! 


= (Sn-lS* + 


■y)h n 



(4.30) 

(4.31) 
(4.32) 

(4.33) 



(4.34) 

(4.35) 
(4.36) 
(4.37) 
(4.38) 
(4.39) 



(see also Appendix A). This 
I) and evolutions ([?.19|), (^.21) 



From the orthogonal polynomials, one obtains the same equations with a = /3 = 7 = 0[ 
differ ence , however, does not look crucial by the follo wing r eason . One can use identification ( 
and (3.23) in order to get the RTC equations (2.19) and (2.20) with u(e) = —1, e = i and 7 = (one can equally compare 
t — i-evolutions)!^]. Let us note that the RTC itself already does not depend on a and fi. Therefore, besides zero 7, one can also 
put a and [3 to be zero. Put it differently, different a, (3 and 7 describes different representations of the same RTC hierarchy. 



5 Forced RTC hierarchy 
5.1 RTC-reduction of 2DTL 

Now we are going to formulate in some invariant terms what reduction of the 2DTL corresponds to the RTC hierarchy. For 
doing this, let us return again to the Lax representation (3.58) embedding the RTC into the 2DTL. Using (3.6), one can easily 



n The standard form of the RTC with an arbitrary value of e and v as in (2.12) can be reached by the proper redefinition of 
time. 



1G 



prove the following identities 

N a c* 

E^fc-l'Sfc-l _ 1 1 
h k h N hjv-i 

k—n 

(5.1) 

N 

y ' SkS k h k = h n — ftjv+i 

fc — n 

Because of these identities, the matrices C and C T have z ero modes ~ S k —i and S k _ 1 /h k respectively. Therefore, one could 
naively expect that they are not invertible and get (using ( |5.l| )) that 

(CC) nk = &nk 



(5-2) 

[££)nk = 4ifc ; 

Since the reduction is to be described as an invariant condition imposed on C and C, these formulas might serve as a starting 
point to describe the reduction of the 2DTL to the RTC hierarchy only if their r.h.s. does not depend on the dynamical 
variables. It seems not to be the case. 

However, these formulas require some careful treatment. Indeed, the formulation of the 2DTL in terms of infinite- 
dimensional matrices, although being correct as a formal construction requires some accuracy if one wants to work with the 
genuine matrices since the products of the infinite-dimensional matrices should be properly defined. In fact, this product exists 
for the "band" matrices, i.e. those with only a finite number of the non-zero diagonals, and in some other more complicated 



cases (of special diver genc y conditions). One can easily see from (3.58) that the RTC Lax operators do not belong to this class. 



Therefore, equations ( |5.2j ) just do not make sense in this case (this is why the interpretation of the general RTC hierarchy in 
invariant (say, Grassmannian) terms is a little bit complicated, see the discussion in the last section). 

Moreover, the presence of the zero mode does not mean automatically that the matrix is not invertible since this zero mode 
is to be normalizable (this is the counterpart of the band structure of matrices, the number of non-zero entries in the column 
describing the zero mode is to be finite, or to satisfy some fast-decreasing conditions).^] 

However, in the case of forced hierarchy, some of the indicated problems are removed since one needs to multiply only 
quarter-infinite matrices and, say, the product CC always exists. Certainly the inverse order of the multipliers is still impossible. 



Therefore, only the first formula in (5.2) becomes well-defined acquiring the form 

(CC) nk = S nk (5.3) 

This formula can be already taken as a definition of the RTC-reduction of the 2DTL in the forced case as it does not depend 
on dynamical variables. Now we will show how this definition is reflected in different formulations of the 2DTL. 



5.2 Fermionic representation 

Here we would like to describe the 2DTL hierarchy in terms of the massless fermions Q and to describe manifestly in the 
subspace in the Grassmannian corresponding to the RTC-reduction. In this and the next subsections we closely follows the 
papers [0 and Q where further technical details can be found. 

Let us define the fermionic operators on sphere 

iP(z) = J2i, k z k , V*(2) = ^Vfc^ fc (5.4) 
feez fcez 
with fermionic modes satisfying the usual anti-commutation relations: 

{tpk, = <5fcm , {i>k, ^m} = {ipt, ip* m } = (5.5) 

The Dirac vacuum |0) is defined by the conditions: 

^, fe |0) = , k < ; Vfc|0> = , k > (5.6) 

We also need to introduce the "shifted" vacua 

i/> k \ri) = , k < n ; xl>l\n) = , k > n (5.7) 

12 It is evident from the following simple example. Let us consider the band matrix B = I — A. This matrix is invertible, but 
its inverse matrix / + A + A 2 + . . . does not possess the band structure. At the same time, there exists the zero mode fj = const 

of the matrix Bij , but it is non- normalizable since the product f T f = ( const ) is divergent. 
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From fcrmionic modes one can built the £/(l)-currcnts 



Jk 



iez 



Mi. 



J-k = Jk , k e Z4 



and define "Hamiltonians" 



H(t h ) = ^t k J k , H(t_ k ) : 



\ " 



t-kJk 



(5.8) 



(5.9) 



fc=i fc=i 

where {t k } and {t- k } are "positive" and "negative" times correspondingly, which generate the evolution of nonlinear system. 
Let g be an arbitrary element of the Clifford group which does not mix the ip- and ip*- modes : 



g=:exp A k 



(5.10) 

where : : denotes the normal ordering with respect to the Dirac vacuum |0). Then it is well known (sec Q and references 
therein) that 

r„(t) = (n\e H ^ge- a ^-^\n) (5.11) 
solves the two-dimensional Toda lattice hierarchy, i.e. is the solution to the whole set of the Hirota bilinear equations. Any 
particular solution dep end s only on the choice of the element g (or, equiva lently , it can be uniquely described by the matrix 
A km ). From relations (Eji) one can conclude that any element of the form (5.10) rotates the fermionic modes as follows 



gipkg 



ipjRjk , 9^* k 9 = ip]R 



(5.12) 



where the matrix Rj k can be expressed through Aj k . We will see below that the general solution (5.11) can be expressed in 
the determinant form with explicit dependence on Rj k . 

Now we deal with the forced hierarchies, i.e. impose the condition 

r„ = , n < (5.13) 
Let us determine what substitutes the general expression (5.1C) in this case. It is reasonable to look for the point of the 
Grassmannian in the form 

g = g P+ (5.14) 

where P+ is the projector onto positive states: 

P + \n) = 9{n)\n) (5.15) 
It can be constructed from the fermionic modes as follows 



=: exp 



(5.16) 



and enjoy the properties 



P+^. k = i>-kP+ = , k > ; [P+,i> k ] = [P+,i>Z] = , k > ; Pl=P+ 



(5.17) 



The insertion of such a projector into eq. (5.11) leaves us with go depending only on ip k and ipf with k > 0. Therefore, it 
has the form [K2l 



90 



: exp ■ 



A(z, u>)i/'-|-(z)i/>i (w 1 )dzdw 



i>0 



(5.18) 



where V+( 2: ) = ^2 k >o^ >kzk ' ^+( z ) = ^2 k >o^k z ~ k anc ^ 7 ' s some integration contour. The matrix elements Aij from (5.10) 
are immediately connected with the modes of the function A(z, w) and, in the forced case, this matrix coincides with the quarter 
of the matrix Rij from (5.12) (the rest part of R g is just unit). Now we are ready to describe manifestly the RTC reduction 
of the forced hierarchy. The reduction condition ( |5.3| ) means that the matrix Aij celebrates the property 

A+-A-A_=A (5.19) 
where A+ and A_ are the quarter-infinite matrices, (A+)jj = <5i,j-i, (A— )jj = 5; 
would be inverse to each other in the case of infinite matrices), although A+ ■ A_ 



.1 j, These matrices are not invertible (they 
= 1. The point is that the inverse order in 



this product does not lead to the unit matrix. By the same reason, condition (5.19) can not be replaced, sa y, b y the similar 
condition A+A = AA+ (or any of remaining two possibilities). This reflects the fact that only the first of eqs.(5.2) makes good 
sense in the forced RTC case. 



Condition (5.19) means that Aij = Ai^j and 



A(z, w) 



2niz 



Siz-w- 1 ) (5.20) 



where £t(z) is some arbitrary measure (compare with ( p.l[ )). 

Thus, RTC-reduction is described by Toeplitz matrices Aij = Ai—j which give rise to the corresponding subspace in the 
whole Grassmannian given by arbitrary matrices Aij. In the next subsection, we demonstrate this by the direct calculation 
from the unitary matrix model representation of the RTC. 
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5.3 Determinant representation 



One can show from formulas ( p,14) )-( p.l8| ) that r-function (5.11) of the forced hierarchy has the following determinant repre- 
sentation |H| 



T n (t) = det 



dl 1 (-d t _ i y j A{z,w)exp\ J2(t 

•'T [m>0 



t Z m - t- m W~ m ) > dzdw 



(5.21) 



— 0, . . . ,71— 1 



Let u s demonstrate now that the unitary matrix model (= forced RTC) leads to this determinant representation with A as in 

(Eaj). 



We begin with rewriting the orthogonality relation (3.2) in the matrix form. For doing this, we define matrices D and D* 



with the matrix elements determined as the coefficients of the polynomials <&n{z) and &„(z) 



(5.22) 



Then, (3.2) looks like 



D-C ■ D* 



H 



(5.23) 

where superscript T means transponed matrix and H denotes the diagonal matrix with the entries Ca = h%—\ and C is the 
moment matrix with the matrix elements 



dfi(z) , 



2niz 



z'-Jexp{ > (t m z m - t- m z~ m ) 



v m>0 



(5.24) 



Let us note that D (D* T ) is the upper (lower) triangle matrix with the units on the diagonal (because of (3.4)). Indeed, this 
represen tation is nothing but the Riemann-Hilbert problem for the forced hierarchy. Now taking the determinant of the both 
sides of (5.23), one gets 



det dj = \ \ h k = T n 

nxn -L J- 

due to formula ( [j.3| ). The remaining last step is to observe that 

d)i(z) 



c« = fi? 1 (-a e _ 1 ) J 



2mz 



exp ■ 



. m>0 



(t m z m ~t„ m z- m ) s =di 1 (-a t _ 1 yc 1 



(5.25) 



(5.26) 



T n (t) = det 



dU-dt^y 



dfi(z) 
2niz 



exp 



E 

.m>0 



{t m Z t — YfiZ ) 



(5.27) 



This expression coincides with (| 5.2l l ) wi th A(z,w) chosen as in ( |5.2C| ). One can also remark that the moment matrix dj is 
Toeplitz matrix satisfying formula ( |5.19| ). This proves from the different approach that the RTC-reduction is defined by the 
Tocplitz matrices. 



6 Relativist ic Toda molecule 



6.1 General properties 



In this section we consider furt her restrictions on the RTC which allows one to consider the both products in (5.2). Namely, in 
addition to the condition ( |xl3| ) picking up forced hierarchy, we impose the following constraint 

T„=0, n > N (6.1) 

for some N. This system should be called N — 1-particle relativistic Toda molecule, by analogy with the non-relativistic case 
and is nothing but RTC-reduction of the two-dimensional Toda molecule [^, [T] 

Now we describe this system in different representations. Let us start with the general two-dimensional Toda molecule and 



discuss which element g in (5.10) describes this hierarchy. Since the Toda molecule is the very particular case of the forced 



13 Sometimes the Toda molecule is called non-periodic Toda ]25|. It is an immediate generalization of the Liouville system. 
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hierarchy, we can look at representation (5.15). Then, the sZ(TV) Toda can be described by the matrix Aij of finite rank TV 
This means that it can be presented as the finite sum 

JV 
k 

where and g^ are arbitrary coefficients. For the kernel A(z,w) ( |5,18| ), this condition is 



(6.2) 



(6.3) 



where f^ k \z) and gW(z) are arbitrary functions. Indeed, these functions (or the sets of coefficients) describe the way how 
sl(N) group can be embedded into gl(oo) (different embeddings are related by the external gl(oo) automorphisms of gl(N)) 
|15[ - this is why the system is sometimes c alled sl(N) Toda molecule. From this description, one can immediately read off the 
corresponding determinant representation ( p\2l| ). 

Indeed, equation (3.28) and condition ( |j,l| ) implies that log to and log rjv satisfy the free wave equation 

dt 1 dt_ 1 log t = d tl dt_ 1 logTjy = (6.4) 

Since the relative normalization of T n 's is not fixed, we are free to choose to = 1. Then, 

T (t) = 1 , r N (t) = x (ti)x(t-i) (6.5) 



where x(*l) an d x(*— l) are arbitrary functions. 2DTL with boundary conditions (3.5) was considered in [|2Sj]. The solution to 
(3.28) in this case is given by [Q: 

r„(i) = det d>- 1 (-d t _ i y- 1 T 1 (t) (6.6) 

with 



n(t) = ^y fc >(t)a w (*-i) 
fc=i 



where functions a( fe '(t) and a( fe '(t_i) satisfy 

dctd'- 1 a^(t)=x(t) , dct(-dt_ i y- 1 a^(t_ 1 ) = x(t-i) 



(6.7) 



(6.8) 



This result coincides with that obtained by substituting into (5.21) the kernel A(z,w) of the form (3.3) 



Let us stress that, although the Toda molecule looks as the forced hierarchy with one mor e pr ojector inserted, this is 
described by the infinite number of the fermionic modes, since infinitely many matrix elements of (6.2) are not zero. However, 
the Toda molecule can be described by the finite matrix Lax representation how we shall see in the next subsection. 

Now we have to describe the RTC-r educt ion of the Toda molecu le. It is clear that, to do this, one needs to impose on 
the Cn-element of the moment matrix (5.26) to have the structure ( |6.6| ), or, which is the same, to require that the matrix 
d li 1 (-dt_ 1 ) j - 1 J2k=i a (fc) (*)o (fc) (*-l) to be Toeplitz. These conditions leads to some equations with the solution just giving 
the RTC Toda molecule, i.e. picking up the RTC-reduction among all the Toda molecule solutions. The same reduction also 
can be formulated as the condition for matrix (6.2) to be Toeplitz which looks quite tricky since naively leads to overfulled 
system of equations. 

However, in the last subsection of this section we investigate the simplest si (2) case in detail and demonstrate that these 
conditions have the non-trivial solutions, i.e. the relativistic Toda molecule does really exist. Indeed, how is already clear 
from the next subsection, the RTC Toda molecule is the system with finite (JV — 1) degrees of freedom, described by TV — 1 
independent time flows and possessing, therefore, finitely many conservation laws. Unlike this, the whole 2DTL hierarchy, and 
its any reductions considered above (RTC, forced etc) are the systems with infinitely many degrees of freedom. The RTC 
molecule is so "small" because of rigidity of the two simultaneous reductions: the RTC and Toda molecule ones. In a sense, 
these two reductions are almost "perpendicular" remaining very small room for common solutions. 



6.2 Lax representation 

In all our previous considerations, we dealt with infinite-dimensional matrices. Let us note that the Toda molecule can be 
effectively treated in terms of TV X TV mat rice s like the forced case could be described by the quarter-infinite matrices. This 
allows one to deal with the both identities (|5.2j) since all the products of finite matrices are well-defined. 



To see this, one can just look at the recurrent relation (3.58) and observe that there exists the finite-dimensional subsystem 



of (TV) polynomials which is decoupled from the whole system. The recurrent relation for these polynomials can be considered as 
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the finite- matri x L ax operator (which still does not depend on the spectral parameter, in contrast to (2.22)). Indeed, from (3.58) 
and condition ( |3,l| ), i.e. h n /h n -i = as n > N (the Toda molecule conditions in terms of S-variables read as S n = S* = 1 for 
n > N — 2 or n > 0), one can see that 



zVn(z) = Vn+i(z) - Vn(z), zVn+i{z) =V N+ 2{z) -V n+ i{z) etc. 



(6.9) 



i.e. all the polynomials V n with n > N are trivially expressed through Vn- Therefore, the system can be effectively described 
by the dynamics of only some first polynomials (i.e. has really finite number of degrees of freedom). Certainly, all the same 
is correct for the star-polynomials V* although, in this case, it would be better to use the original non-singularly normalized 
polynomials 



Now let us look at the corresponding Lax operators ( 4.17 )-( 1.18). They are getting quite trivial everywhere but in the left 
upper corner of the size N X N. For instance, 



So 



1 



— — —Si ~S\Sq 
no 



-1 1 



-1 1 



-1 1 



V 



(6.10) 



,N- 



and analogously for the Lax operator C. Therefore, one can restrict himself to the system of N polynomials Vn , n = 0, 1, . 
1 and the finite matrix Lax operators (of the size N X N). 

Now one needs only to check that this finite system still has the same evolution equations J4.2jj ) . It turns out to be the 
case only for the first N — 1 times. This is not so surprising, since, in the finite system with N — 1 degrees of freedom, only 
first N — 1 time flows are independent. Therefore, if looking at the finite matrix Lax operators, one gets the dependent higher 
flows. On the other hand, if one embeds this finite system into the infinite 2DTL, one observes that the higher flows can be no 
longer described inside this finite system. Let us remark that just the finite system is often called relativistic Toda molecule 
(see, for example, Jisj ) ■ ) 

To simplify further notations, we introduce, instead of S n , S*, t he new variables s n = (— ) n+1 S n , s„ = (— ) n + 1 S*. Then, 



one can realize a very interesting property of the Lax operator (3.10) - it can be constructed as the product of simpler ones: 

C = C N C N - 1 ...C 1 (6.11) 
where Ck is the unit matrix wherever but a 2 X 2-block: 



/ 



Analogously 
with 



\ 



i J 

C = Ci . . . £jv-i£jv 
\ 

G k = ( 

1 J 



s k s k - 1 



1 - SUS 



(6.12) 



(6.13) 



(6.14) 



One can trivially see that C k C k = C k C k = 1, and, therefore, from (3.11) and (3.13), one obtains ££ = ££ = 1 (cf. fQ) 
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From formulas ( 3.11 )-( 5. 13 ) , one trivially obtains that det C = det C = 1 which reminds once more of the sl(N) algebra. 
More generally, the factorization property of the Lax operators opens the wide road to the group theory interpretation of the 
RTC molecule. Indeed, following the line of papers [ p^ , [l5| |l(| , one should identify the (family of) solutions to the integrable 
hierarchies with representations of the algebra of functions on the underlying group manifold. In the quantum group case, one 
should expect that different reductions correspond to fixing the irreducible representations |l4[ Fil . As the classical counterpart 
of this statement, one should consider the (quadratic) Poisson structure for the group elements 14 

{g ® g} = [ns® a] ( 6 .i 5 ) 

where r is the classical r-matrix (see, e.g., ||), and, instead of irreducible representations, look at the simplectic leaves, i.e. the 
submanifolds where Poisson structure is non-degenerate (an analog of the geometric quantization) - see Jlt| ). This procedure 
has been performed in ]l6[ for the sl(N) group and was shown to lead to the relativistic Toda molecule for the leaves of the 
dimension equal to the doubled rank of the group (i.e. 2(jV — 1) for sl(N)). 

It was shown in |l^ that there always exist such leaves^] and they were manifestly described. Namely, for the sl(2) case, 
the 2-dimensional simplectic leaf can be described by the group element of the form 

ad - b 2 = 1 (6.16) 




Multiplying this matrix by the trivial rescaling matrix g = ^ ^ J one can easily transform g to the matrix Gk in ( |6.12| ). 

In fact, the only important information of this matrix at the moment is that it is parametrized by a 2-dimensional manifold. 
However, the crucial test comes when considering the structure of this matrix for higher rank groups. Namely, for arbitrary 
sl(N) group, the group element corresponding to the proper simplectic leaf is constructed from elementary sl(2) building blocks 
exactly how it is done in formula (6.11) Jl7|, |l(j]. Therefore, the Lax operator of the relativistic sl(N) Toda molecule can be 



treated as the group element of the special form corresponding to the simplectic leaf of the dimension 2(N — 1). Then, one can 



construct the Hamiltonians commuting with respect to the Poisson structure (6.1!:) and giving the time flows just as traces of 
the Lax operator. 



It may seem that there does not remain any room for the other Lax operator (3.13). However, this is not the case. Indeed, 
the above construction is correct up to any Weyl transformation. The transformation corresponding to the longest element 
of the Weyl group maps C operator to C operator. This explains the symmetry of dynamics with respect to the positive and 
negative times. As for other Weyl transformations, the corresponding Hamiltonians can be hardly into 2DTL dynamics. 



6.3 sl{2) molecule 

To illustrate the results of this section, let us consider the simplest sl(2) example of the relativistic Toda molecule. This case is 
to be considered as a "relativization" of the Liouville theory. Surprisingly enough, this case turns out to be equivalent to the 
Liouville theory. This is a specific feature of the sl(2) case, which is no longer correct already for the sZ(3) group. This can 
be demonstrated by comparing the Lax operators, but we use more immediate way of comparing solutions to the equations of 
hierarchy. 

In fact, the sl(2) hierarchy corresponding to the rank one group contains only two equations: one for positive and the other 



one for negative times. One of these equations is (1.1) with the Toda molecule constraints g_i — * oo, q2 — + -co etc. Under 
these constraints, there are only two nontrivial RTC equations 

e <31-<30 



ijo = (1 + e?o)(l + eqi) 



91 = -(1 + e?o)(l + eqi)- 



1 + e 2 eii~io 

91-10 



(6.17) 



1 + e 2 eii-io 

These equations imply that go + gi = C\t + C2 with arbitrary constants Ci, Ci- Now introducing the new variable q 
go — ^(C\t + C2) one gets the equation 



(1 + \e Cl ) 2 - e 2 q 2 



-2q 



(6.18) 



l + t 2 e- 2 i 

14 In paper ]26[, the same ideology was applied to the linear Poisson structures. This leads to the non-relativistic TC, in 
accordance with the "algebraic" character of this latter as opposed to the "group" character of the RTC. 
15 When quantizing, they lead to the infinite-dimensional representations of the algebra of functions. 
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The limit of e — > leads to the Liouville equation 



-2q 



Equation (3.18) has the solution 



dq 



(l+ieOi) 2 +0 3 



ii +C 4 



■C 3 e- 2 i 



(6.19) 



(6.20) 



(with some redefined constants {C;} 
(non-relativistic Toda molecule) 



i.e. q = log [A Q sinh(Ai*i + A 2 )\ 
{Ai\) which is to be compared with the analogous solution to the Liouville equation 



dq 



h + C 4 



(6.21) 



V C 3 ~ e~ 2q 

One easily observes that these both can be transformed to each other by the trivial rescaling of the time variable and integration 
constants. Therefore, "the relativization" of the Liouville equation does not lead to a new equation. 

Let us not e, h owever, that switching on the negative time changes this statement. Indeed, looking at equation completely 
equivalent to (1.1), one can check that the solutions really depend on an arbitrary linear combination of times: 



log [Aq sinh(Aiti 



■A 3 t- 



(6.22) 



Therefore, the notorious property of dependence only on the sum of positive and negative times, which distinguishes the TC 
hierarchy is not fulfilled for the RTC molecules. Thus, these systems are really different even in the sl(2) case. 



In fact, the pair of equations (6.17) can be really substituted by the only equation for the r-function t\ (3.28) with the Toda 
molecule ( |S.^ ) and reduction conditions. To do this, let us first note that conditions (6.4) allows one to normalize r- functions so 
that to = Tjv = 1. In particular, in the sl(2) case this means that one can add to the equation (3.28) for ri analogous equations 
for to and T2 with the coefficient to get 



-dt 1 d t _ 1 log 

2 T T 2 



TQT2 



and, for the new function r = 



/TQT2 ' 



dt 1 d t _ 1 log 



Using formulas (3.9) and (3.13), one obtains that 



-2q 



(6.23) 



(6.24) 



(6.25) 



We have already remarked (see the very end of sect. 3.1) that, generically, the connection of r-functions with g-variables may 
require some additional shift of these latter. 



Now one should solve equation (3.24) with the moment matrix Q\ 1 (— 9t_ 1 ) J 1 Ti(t) constrained to be Toe plitz. In our 
case, this means that dt 1 dt_ 1 T = — t. Then, one obtains the equation (the same equation can be obtained from ( p. 18| ) by the 
substitute (3.25)) 

r 2 - 1 



(M 



This equation has the solution 



(h 



B 2 



(6.26) 



(6.27) 



t = sinh(Biii + B 2 ) 



(6.28) 

and can be reduced to (6.20) through the replace (3.25) r = i e 9+const This formula differs from (3.25) by a constant shift of q. 
Indeed, we have already remarked (see the very end of sect. 3.1) that, generically, the connection of r-functions with q-variablcs 
may require some addition al sh ift of these latter. In the present case, the shift depends on value of C'i and is in charge of the 
non-unit coefficient Aq in (6.22). The correct connection of r and q restores if Ci is chosen to be — — , i.e. q n — + q n — -t\. 

Now let us look at the dependence of r on the negative time. Naively, using ( |6.28[ ) and constraint dt 1 dt_ 1 r = — t, one gets 
that the general solution is given by the formula 



: sinh(i?iti 



B 2 ) 



(6.29) 



which contradict s to ( 3.22 ) with arbitrary unrelated coefficients in front of times. This puzzle is solved by noting tha t the 
moment matrix (5.26) can be multiplied by arbitrary triangle matrix with units on the diagonal so that the determinant (3.27) 
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does not change. This means that one really needs to consider the weaker Toeplitz condition dt 1 dt_ 1 T + adt_ 1 T - 
a is a constant (this does not effect equation (5.26)). This results to the general solution for r-function 

r = sinh^iti + B.T 1 *-! + B 2 ) 



where 



(6.30) 



consistent with (3.22). This provides us with some nontrivial example of the simultaneous solutions to the constraints of the 



relativistic Toda molecule and Toeplitz moment matrix. The measure //(z) ( p.20| ) can be obtained from the explicit formulas 
i.2q) and (5.27) by the Fourier transform. 



7 Discrete evolutions and limit to Toda chain 
7.1 Darboux-Backlund transformations 

In this section, we are going to discuss some discrete evolutions of the RTC given by the Darboux-Backlund transformations 
and their limit to the usual Toda chain. One can easily take the continuum limit of the formulas of this section to reproduce 
the TC as the limit of the RTC, both with the standard continuous evolutions. 

The discrete evolution equations in the RTC framework were recently introduced by jj| in a little bit sophisticated way. 
Here we outline the simple approach based on the notion of the Darboux-Backlund transformation (DBT). More details will be 
presented in the separate publication Jlsj ] . 

Let discrete index i denote the successive DBT's. The spectral problem now can be written as follows: 

$ n+1 (i\z) + a n {i)<& n (i\z) = z j$ n (i|z) + b„(i)<I> n _i (i\z) j (7.1) 

In fact, we consider two pairs of different DBT's, the two in each pair being complimentary to each other. This is why we call 
the DBT's in a pair forward and backward DBT's. 

Let us define the first forward DBT (treating it as the discrete evolution) in the form very similar to that of the usual Toda 
chain: 

0„(i + l|z) = <Mi|z) + o4 1) (i)*n-i(i|z) ( 7 - 2 ) 
where QL„ (i) are some unknown functions. From the point of view of the whole hierarchy, this new evolution means nothing 
but a discretized version of the first negative flow equation ( [2.2[ ). One requires that "3? n (i + 1) satisfies the same spectral 
problem as ( |7.l| ) but with the shifted value of i : 

*„+i(i + l\z) + a n (i + l)$n(i + l|z) = z j$„(j + l\z) + b n {i + l)$„_i(i + l\z)\ (7.3) 
Then the compatibility condition gives the equations of the discrete RTC: 



a n (i) — ail, (i) 

a n (i + l) = a n -i(i) — ^y- (7-4) 



On-lCO - «n (0 

b n (i + l) = b n - 1 (i) MO-o^ffl 

Zi ~7TT~ = Zi ~ a "W + a L+i(*) (7.6) 
aj, ' (i) 

where Zi are arbitrary constants. These equations give some natural generalization of the corresponding discrete evolution 
for the Toda chain hierarchy; moreover, there exists the simple limit to the usual discrete Toda equations (see the discussion 
of this limit below). Using representation ( |2.6[ ), ( |2,7| ) it is very easy to rewrite the system ( |7.4| )-( |7^ ) in terms of "coordinates" 
q n {i) only [|. 

From the continuum picture, one can get that there is another simple evolution (ti-flow) treated in a parallel way with 



that generated by the t_i-flow. Therefore, one can try to find the corresponding discrete analog of (2.16). Indeed, such analog 
exists and has the natural form 

* n (i + X\z) = (l - o4 2) (i))<Mi|z) + 2 o4 2) (i)<I> n _i(i| Z ) (7.7) 

(2) 

where Oi n '(i) are some new unknow n fu ncti ons of the corresponding discrete indices. We refer to this evolution as to the 
second forward DBT. Substitution of ( |r.7| ) to (7T) gives quite different system of the discrete evolution equations: 

1 - Cc' 2 ' (i) 

a n (i + l) = a n (i) ^±i^ ) (7.8) 
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b n (i + 1) = &n-i(i) 



1-Cti 2) (r) 

ai 2) (i) 



(7.9) 



(7.10) 



This syste m is the discrete counterpart of the continuum system (2.17), ( |2,18| ). Written in the terms of q n (i) , equations 
( 7.8 )-( 7. 1C ) define the same evolution of coordinates as in the case of the first forward DBT. 

Actually, in four different discrete systems of the RTC equations were written. From our point of view, the additional 
evolutional systems result from the backward Darboux-Backlund transformations which are complimentary to those 
described above. For example, the first backward DBT has the form 



* n (i|*) = -J— { (l - /3« (i)) <fr n+1 (i +l\z) + zP™ (i)<£n(i + l\z)\ 



(7.11) 



where z; are the same constants as in (7.6) and variables /3„ (i) are related with the variables of the first forward DBT as 
follows: 



P^ii) 



(7.12) 



The system obtained from (7.11) can be easily reduced to ( |7.4[ )-( |7.6| ) and is not considered here anymore. 



7.2 Continuum limit 



Introducing some discrete shift of time A > , one can rewrite all the equations describing the first forward DBT as follows: 



*n(t+ A) =* n (f) + o£ 1) (t)# n _i(*) 

a n (t + A) = a n -i(t) 



fe n (* + A) = 6„_l(t) 



a„_i(t) - a^t) 
&„(*) -a^W 
^-iW-ai'l,!*) 



(7.13) 
(7.14) 

(7.15) 



In order to get the proper limit in (7.6) one should rescale the constants Zj ; for simplicity, we assume now that they do not 
depend on i : 

Zi — » g A 

Then we get 



and solution to (7.17) has the asymptotics 



a«(t) = - ff A^# + 0(A 2 ) 



a„(t) 



and, therefore, from (7.13)-(7.15) 



ij$„0,t + A)-$»(z,t)j 
i(a n (t + A)-a„(t)l = - 9 f-^ 



9 ^ 4n _ l(2jt ) + (A) 
a n (t) 

&n(t) 6n+l(<) 



O(A) 
- O(A) 



a n+ i(t), 

l{5 n (i + A)-5 n (t)) = - fl 6„(i)f 3_ y 

and, in the continuum limit, A — » these equations lead exactly to (2.4), ( p.q ). 

The analogous equations can be written for the second forward DBT but now Zi should be rescaled as follows: 

1 

z i ~> 7" 

9 A 



and (Ej), - (7. 10) now have the form 

*„(t + A) = (l-al 2) (i))<I- n (t) + zQ:l 2) (i)<I>„-i(t) 



Substituting to (7.1) gives the equations: 



1 - a <2) (t) 

a n (t + A) =On(t) ^±f— ) 

l-ai 2) (t) 



(7.16) 
(7.17) 



(7.18) 

(7.19) 

(7.20) 
(7.21) 

(7.22) 
(7.23) 
(7.24) 
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b n (t + A) = 6 n _i(t) 



a& 2) (t) 



an(t) + 6„(t) 



l-al 2) (t) _ 1 



From (7.2£), it follows that 



and, therefore, 



6n (*) 1 1 + 9 A f a„ (t) - b„ (t) - b n+1 (t) J | 

ij*„(jM + A)-* n (*,t)} = -SM*)(*n(2,*)-2*»-l(2,t)j + O(A) 
i|on(t + A)-o„(t)| = 9d„(()|l)„(()-ii„ + i(()J + O(A) 
-{& n (tH-A)-&»(t)J = 5b n (t)^a„(t)-o„_ 1 (t) + 6 n _i(t)-6„ + i(t)J + O(A) 



It is clear that, in the limit A — > , one reproduces the continuum equations (2.17), (p.lq) 



(7.25) 
(7.26) 

(7.27) 

(7.28) 
(7.29) 
(7.30) 



7.3 Limit to Toda chain 

Here we very briefly outline the theory of the discrete evolution of the usual Toda chain (see, for example, [ gTJ and references 
therein). The spectral problem for the Toda chain has the form 



Atf„(i|A) = * n +i(»|A) - p„(i)*n(»|A) + ^(i^n-iM*) 



(7.31) 



where we consider the whole set of successive DBT labelled by the discrete index i as a discretized variant of the continuous 
evolution equations 

5*„(i) = #„(i + l) - * n (t) = A n (i)*„_i(i) (7.32) 



which should be compatible with the spectral problem (7.31). (7.32) describes the forward DBT. Analogously, the backward 
DBT is 

*„(i|A) = ^-^-|* n+1 (i+l|A)+Bn(i)*n(i + l|A)| (7.33) 

where coefficients B n (i) should be determined. The compatibility conditions give the following system of the evolution 
equations: 

p n (i) = -\i - A n+1 (i) - B n (i) (7.34) 

p n (i + 1) = -Ai - A„(i) - B n (i) (7.35) 

R„(i) = A n (i)B n (i) (7.36) 

H n (i + l)=A n (i)B n _ 1 (i) (7.37) 

where A; serve as free parameters. In order to get the close connection with the general theory of integrable systems as well 
as with the theory of orthogonal polynomials, one can introduce new variables h n (i) = exp(q n (i)) through the relation: 

h n (i) 



Then, in the terms of variables h n (i) 



Rn(i) 



A»(i) 



B n (i) 



hn-l(i) 

h„(i + l) 
h„-i(i) 
h n {i) 



(7.38) 



(7.39) 



(7.40) 



M« + l) 

and, from (7.34)-( fr.35| ), one can get the following evolution equations for the discrete Toda chain (assuming for simplicity that 
parameters A^ do no depend on i): 

h„(i + l) 



hn(i + l)hn(i — 1) hn-i(i — 1) 



hl(i) 



1 



h n+1 (i + 1) 
h„(i - 1) 



(7.41) 



Now we can consider the transition from the RTC to the usual discrete-time Toda chain. Let (compare with (2.6), (2.7)) 

<Xn(i) ~ 1 — ep n (i) ; b„(i) ~ -e 2 R„(i) 
z~l + eA; Zi~l + eA; 



(7.42) 
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Introduce also functions 4f n (i)r 



(7.43) 



(7.44) 



It is easy to see that (7.2) leads to the forward DBT for the Toda chain if one identifies 

a£\i) ~ eA n {i) 

In this case, for example, equation (7.6) reduces to (7.34). Actually, it is easy to see that the whole system j7.4|) -( p!q ) reduces 
to (7.34)-(7.37) in the limit e — > . Moreover, from (7.12) one gets 



R n (i) 
' A n (i) 



+ 0(e 2 ) ~eB n (i) 



(7.45) 



thus obtaining the proper limit of the first backward DBT of the RTC (7.11) to that of the Toda chain (7.33). 

It is clear also that, in th is limit , the second forward DBT ( \i.t\ ) reduces to ( 7.32 ) and all the equations we get from (7i 

(7.10) are equivalent to ( 7.34 )-( 7.37) when e -> 0. 



There exist some other interesting limits. In particular, one can consider the degenerate case of evolution equations (7.4)- 
(7.6) in the limit z; — > oo. In this case, the solution to (7.6) has the asymptotics 



Thus, in this limit equations (j^4|) and ( |7.5| ) take the simple form: 



a n {l + 1) = On-x(t) 7TT — TTT 

On-l(l) - b re (l) 

6n+l(i) - 



6n(« + 1) = b n (i) 



b„(i) - a„-i(i) 



(7.46) 



(7.47) 



(7.48) 



We should stress that these equations can be easily derived from the spectral problem (7.1) just representing it in the form 

<£„ +1 (i|z) + a„(i)$„(i|z) = z $ n (i + l|z) (7.49) 
*„(i + l|z) = *„(i|z) + fe„(i)<I> n _i(i|z) (7.50) 



Using representation ( |2.6[ ), ( |2 . 7| ) we can rewrite system ( |7.47| ), ( |7.48[ ) in terms of coordinates q n (i) as follows: 

1 + e 2 exp{q n+ i(i) - q n (i)} 



cxp 



{q n (i + 1) - 2q n (i) + q n (i — 1)} 



1 + e 2 exp{g„(i) - g n _i(i)} 



(7.51) 



thus getting the counterpart of the discrete-time Toda chai n j3[ . We should stress, however, that equations (7.41) and (7.51) arc 
obtained as different limits from the general discrete RTC ( |7.4| )-( p^ ). Moreover, they are described by different Lax operators. 



8 Some applications to biorthogonal polynomials 

In this section we discuss some peculiar properties of the system of (rclativistic) polynomials leading to the RTC, biorthogonal 
with some specific measures. 



8.1 Finite systems of orthogonal polynomials 



(8.1) 



We already discussed in sect. 3 that the relativistic polynomials satisfy recurrent relations (3.5) 

#„+l(z) = Z® n (z) + SnZ" 1 ®^- 1 ) 

provided the polynomials are normalized to be the n-th order monic (i.e. with the coefficient in front of the leading term equal 
to 1) polynomials. In particular, the initial conditions are $o(^) = &*)( z ) = 1. We also demonstrated that such polynomials 
satisfy the RTC recurrent relations (2.1) 



#„ + i(z) + a„$ n (z) = z{*„(z) + 6 n 4„_i(z)} , n e Z 



(8.2) 

16 Let us note that, under this limit, the cor responding biorthogonal polynomials turn into the usual Toda orthogonal poly- 
nomials, with the biorthogonality relation (3.2) transforming to the orthogonality relation. This is why we sometimes call these 
biorthogonal polynomials relativistic ones. 
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As far as we know, it was G.Baxter who first considered the system (p.5[ ) [|28j (for recent development of the theory of these 
polynomials see e.g. |29[ On the other hand, Pastro showed |29| that the theory of the polynomials defined by (3.5) is 



closely related to the Laurent orthogonal polynomials first introduced in the framework of the strong Stiltjes moment problem 
[32) . It is interesting to note that the RTC was studied in J33) by the method of continued T-fractions arising in the context of 
the same strong Stiltjes moment problem, however, without using the orthogonal polynomials. 

Note that in general, S n and S* are arbitrary complex parameters. When S n = S* and \S n \ < 1, one gets the theory of 
Szego polynomials orthogonal on the unit circle ^5[ . The case of \S n \ > 1 is much less trivial and has not been considered 
in detail. 



It is interesting to discuss how given recurrent relations ((3.5) or (p.l[)) are connected to the concrete orthogonality measure. 



In general, when the space of the polynomials is infinite-dimensional, one can hardly formulate any statement on the existence 
and support of the corresponding measure. However, for the finite-dimensional representations corresponding to the Toda 
molecule such (discrete) measures can be found explicitly. 



For doing this, note that one can easily derive from (3.5) the identity 



»-»(*;(»)*„+!(*) - **(*)*„+i (vm- 1 = (x -y)^2 (y -1 )'* 1 (8.3) 

where ho = 1 and we introduced for the sake of brevity $£(2) = z n <^^ l (l/ z). This identity is valid for arbitrary complex x and 
y unless x ^ y. For x = y we have the identity 

71 

(y n h n r 1 @n(,y)K+i(y) - (2/)) = ^WKiVvfri 1 (8.4) 

k=0 

Identities ( BJ3 ) and can be considered as "relativistic" analogs of the well known Darboux-Christoffel formula in the theory 
of the ordinary orthogonal polynomials jji) . 

Now let us turn to the finite set of polynomials. In order to reduce our system, we put h n +l =0, h n ^ (this means 
that SnS* = I and this corresponds to the Toda molecule case of section 6). Let j/j, j = 0, 1, , n be the roots of the 



polynomial $ n +i(y). We assume additionally that all these roots are simple. Then, from (S.3) and (8.4), one gets the following 
orthogonality relation 

n 

Y&kiVs^UVr 1 )^ 1 ='W- 1 S rs (8.5) 
fc=0 

where 

WS = Wfew (8 ' 6) 

At last, using that, for the orthogonal matrices, their transponed are orthogonal too, one finally obtains the dual orthogo- 
nality relation 



Wk<&n{ykW m {y k 1 ) = h n Snm (8.7) 



fc=0 



Therefore, uij. is the discrete weight function whose support is the set of zeros of the polynomial i(z). 

Note that because of the condition S n S* = 1, the roots y^ of the polynomial ( I > n +i(2) coincide with the roots of the 
polynomial < I > *^ 1 (l/z). Then, using the symmetry between the polynomials Q n (z) and <E>* (2), one derives another orthogonality 
relation 

n 

'wfc'Mj/^^fefe) = h n 8 nm (8.8) 



fc=0 



where 



™k = - . rz~, - —. (8.9) 



In general, the weight functions w k and Wj, do not coincide and, therefore, we have two different orthogonality relations for 
the same polynomials. However, if Sk = k = 0, 1, . . . , n (i.e. we deal with the symmetric case) then = w^. Moreover, 

in this case the polynomials $fe(z) and &^(z) coincide and, hence, the set of the roots y k of the polynomial "Iv+i coincides with 
the set y^ 1 of roots of the polynomial < I > *^ 1 (z). This means that, in the symmetric case, the polynomial <I> n _|_i(2) is invertiblc, 
i.e. the coefficients in front of the terms z-> and 2™ +1 ~ J coincide (or, equivalently, all the roots enter with the pairs y^, V^ 1 )- 
In particular, when \Sk\ < 1, all the roots lie on the unit circle. 



28 



8.2 Anzatz of separated variables 

Now let us consider another interesting peculiar case of the relativistic polynomials connected with a specific solution to the 
RTC hierarchy. Namely, look at the following ansatz of separated variables (see (2.1)) (for the non-relativistic Toda chain this 
ansatz was studied in 

On(t) = 7l(*K, b n {t) = 72(i)Mn (8.10) 
where fi(t) depend on t only whereas u n and fj, n depend only on n. Putting bo = (because of one gets the solution 

t , - tn 

a„(t) = -, b n (t) = -- — (8.11) 

n + c+1 (n + c)(n + c+l) v ; 

where a is an arbitrary constant. Then, the corresponding coefficients S n and S* are equal to 



a„(t) = — - 



(c + l)„+i 

where (c) n = c(c + 1) . . . (c + n — 1) is the Pochhammer symbol. 



p n {t) = -t- n {c) n+1 



(8.12) 



In order to find explicit expressions for the polynomials <& n {z) and <J>*(z), note that the recurrent coefficients (3.12) or 
(|8,ll| ) can be obtained as some limiting case of the system of biorthogonal polynomials on the unit circle proposed by R.Askey 
|18| . Indeed, the Askey polynomials have the recurrent coefficients 



( b )n + l 
(c+ l) n+ l ' 



(c)n+l 

"(6 + 1) 



n + 1 



(8.13) 



where a, b are arbitrary parameters. Corresponding polynomials are expressed in terms of the Gauss hypergeometric function 

Hi 



(&)> 



(C + l)r, 



— n, 1 + c 
1-n-b' 



(C)r 



(b + l)r 



f-n,l 
VI - n 



+ b . 

— n — c 



These polynomials are biorthogonal on the unit circle 

<S> n (e W )^* m (e- ie )e- i6 ^- b '>/ 2 \ sin(6/2)\ c+b de = 0, n^m 



(8.14) 



(8.15) 



It is easy to see that the recurrent coefficients (S.12) can be obtained from (S.13) by the trivial rescaling in the limit b — > oo. 
Omitting the simple technical details, we present here only the result. The polynomials corresponding to solution (3.11) or 
( J8.12j ) are 



(c+l)„ 

These polynomials are biorthogonal on the unit circle 



2 F {-n,l + c;-z/t), $* (z; t) = t~ n (c) n 1F1 ( " ;tz) 

V 1 — n — c ) 



*"(e i6 W m (eT iB )e ic6 exp(-t e - ie ) = 0, n £ m. 



(8.16) 



(8.17) 



Now let us consider the first Darboux-Backlund transformation (7.47)-(7.4£ ) (see the previous section). We define w(z) to 
be a weight function for the biorthogonal polynomials 

(8.18) 



<$>n{z)<$>* m {z 1 )w(z)dz = 0, nj^m 

Then, it is transformed under the first Darboux-Backlund transformation as 

w(z; t + 1) = kzw(z; t) 

where k is a normalization constant. 



(8.19) 



In particular, in the case of the separate variables solution (8.11), this transformation is equivalent to shifting the parameter 

n the bo 

a n (t; r) 



c: c = c + 1, and, thus, one has a whole family of solutions depending on the both continuous and discrete times 

tn , t 



b„(t;r) 



(n + co + r)(n + cq +t + 1) 



n + co + r + 1 



(8.20) 



The property (S.19) is obviously fulfilled, which can be seen from orthogonality relation (8.17). 

It is interesting to note that the Askey polynomials (8.14) themselves obey the discrete time Toda dynamics determined by 
the Darboux-Baklund transformation provided that c(r) = co + r, b(r) = bo — r. 

The general Darboux-Baklund transformations of the biorthogonal polynomials will be considered in a separate publication. 
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9 Concluding remarks 



In the present paper we considered some different representations of the relativistic Toda hierarchy, which are naively non-related 
to each other. However, from the point of view of studying the RTC hierarchy itself, the most promising representation is that 
describing the relativistic Toda hierarchy as a particular reduction of the two-dimensional Toda lattice hierarchy. However, 
even this quite large enveloping hierarchy is still insufficient. 

Indeed, the RTC-reduction of the 2DTL being transparent on the level of the Lax operators can not be explicitly written 
in terms of the Clifford element. In other words, no reasonable description of the point of the Grassmannian is possible in the 
context of the usual (one-component) Toda lattice. Loosely speaking, the Toda lattice is too "rigid" to reproduce both the 
continuous and discrete flows of the RTC. Therefore, one can try to embed the RTC in a more general system which admits 
more natural reductions. 

This is done in the forthcoming publication Jl9t , where we show that the RTC has a nice interpretation if considering it as a 
simple reduction of the two-component KP (Toda) hierarchy. In this case, any dynamical variables acquire two discrete indices 
(we have two fermionic vacua in this picture) so that one can treat the second vacuum number as a discrete time which generates 
the sequence of the degenerate Backlund transformations (or, equivalently, the discrete evolution of the modified Toda-type). 
On the other hand, we need yet another discrete index in order to generate the whole discrete-time RTC. It turns out that 
the natural interpretation of this additional discretization is possible if one introduces the so-called Miwa variables (with the 
corresponding multiplicities) after imposing the reduction conditions. The whole evolution of the RTC can be interpreted 
now as follows: the evolution along the additional vacuum generates the modified Toda equations while the evolution with 
respect to the multiplicities of the Miwa variables leads to general Darboux-Backlund transformations. This is the reason why 
the (properly) reduced two-component KP (Toda) is the true framework for description of RTC. Moreover, it turns out that, 
in the framework of the 2-component hierarchy, the continuous AKNS system, Toda chain hierarchy and the discrete AKNS 
(which is equivalent to the RTC how we proved in this paper) can be treated on equal footing. In particular, this means that 
all these systems corresponds to the same subspace in the Grassmannian. 
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Appendix A. Evolution equations from the orthogonality condi- 
tions 



Here we shall give some examples of deriving the evolution equations in the polynomial case starting from the orthogonality 
condition ( |j.l| ), ( |3.2| ). In order to keep the close connection with the general theory of the monic orthogonal polynomials, we 
write down here the spectral problem directly in terms of <E> n and <£>*: 



S* 

z<S> n (z) = $ n+1 {z) - Snhn V -^i$ fc (z) = C nk $ k (z) 
fc=0 



(A.l) 



iT 1 **^" 1 ) = ** + 1 ( Z ~ 1 ) - S*h, 



* — ' h k n k 

k=0 



Now the differentiation of (3.2) with respect to t_i gives (for n > k) 



<9i-i ' 



Using (A.2) it is easy to see that 



d$ n hn * a„ 

tt, — = 7 *n-i = (1 - SVi-iS£_i)$™-i = 7-*™-i 

at-i h„-i o n 



(A.2) 



(A.3) 



(A.4) 



30 



Similarly, the differentiation of ( [j,2| ) with respect to ti gives (for n > k) due to (A.l) 



Therefore, 



\ dti ' k 

9*n 



< Z$„,$* k > = SnS^hn 



n-1 ^ A 

fe=0 



(A.5) 



(A.6) 



Sn — 1 hti — 1 

The equations (A. 4), (A.6) are nothing but the simplest evolution equations ( ^.2| ), (2.16) which determine the RTC. Note that 
all the general equations (1.28) can be obtained in the same way. For example, the differentiation of (3.2) for n > k with 
respect to t m gives: 

iw 1 ' ^) = ~< (£ m U<£s, #J >= -(C m ) nk h k e(n — k — 1) (A.7) 

and, therefore, 

d*n ™ , , (A8) 



Similarly, the differentiation of ( [j.2| ) for k = n yields 



(A.9) 



In the particular case of m = 1 , one gets just Cn 



-S n S*_ 1 , which leads to (3.23) with 7 = 0. 
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Appendix B. Lax operators in 2DTL framework 



In this appendix we are going to represent the spectral problem (3.7) arising from the orthogonal polynomials as the "usual" 
spectral problem. First of all, in the fast-decreasing case 



Sn -> , S* -> ; n -> ±oo 
one should prove that the solution to (3.7) with asymptotics (4.19) satisfies the spectral problem j^j 

C-nkVk(z) = zV n (z) 



where C nk is given by (4.17). Indeed, the equation (3.7) can be re-written in the form 



zV k {z) 



— h — (r k (z) 



- zVk-i(z) 



-Vk{z) 



(B.l) 
(B.2) 

(B.3) 



In order to get (B.2), one should sum (B.3) from k = -co to k = n provided the following boundary condition is fulfilled (using 
that h k —* Const when k — > — oo) : 



(k(z) = ^-(v k+1 (z)-zV k (z)j ^0 ; fc^-oo 
With the help of ( |4.2l| ), (juj), it is easy to see that coefficients Wi(k) satisfy the recurrent relations 

— |wf + i(fc + 1) - tUj + i(fc)| - - 1 ~ w i( k ~ !)} = Sl_ 1 w i (k - 1) ; w (k) = 1 



(B.4) 



(B.5) 



^-|u)i(fc + l) -Wj(fc)} 



>Jgfc-jWi-j(fc- j) 

3=1 



(B.6) 



Thus, for any /ired i , the r.h.s of (B.6) tends to zero provided S£u>j(fc) — > as fc — > — oo. For sufficiently fast-decreasing Sk, 

S£ , it is really true. Therefore, we see that all the coefficients of Ck( z ) in (B.4) vanish in the limit k — » — oo. Since solution 

V^\z) is defined in the range of large value s of the spectral parameter z , the whole function Cfe( 2 ) is conv ergent in the limit 

of large negative k and, thus, satisfies (B.4). Note that the same arguments fail in the case of solution ( |4.2C]| ) since it is defined 

at small values of the spectral parameter and the corresponding function Ck( z ) can not be properly defined as k — > — oo and, 

(21 I J 

therefore, V n (z) does not satisfy (B.i) in general. 

To understand this situation in more details, it is instructive to consider the case of the forced hierarchy when 

S k = = 1 ; k < 

S k , Si -> ; k — > +oo 



(B.7) 



The first condition means also (see (3.6)) that 



h 



k + l 



k < h ^ 



(B.8) 



In the forced case, Vn\ z ) can be easily described. Indeed, it is a polynomial solution to (3.7) of the type (3.4) at n > , 
Vq{z) = 1 and 

n < (B.9) 



V£\z) = (z + IT 



and it is well defined in the whole region of the spectral parameter. This solution is non-singular in the sense that 

-^-V^iz) = , n < 



Thus, from (B.3) in the forced case, one can get the same spectral equation as in the fast-decreasing case, i.e 

St. 



zV n (z) =V n+1 (z) - S»h n V -^-V k (z) = C nk V k (z) 

h k 



(B.10) 



(B.H) 



where the summation actually goes from to n due to ( |B.lo| ). 

On the other hand, the second independent solution to ([^t]), Vn ' 1 (z) , is non-polynomial even in the forced case and can 
be treated as a formal series at small values of z only. Due to asymptotics (4.2C), it is singular, i.e. 



-LvL 2) (z)^o n<0 



(B.12) 



17 See identification (3.56) 
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Moreover, putting 



one can find out the exact answer 



Vi 2) (z) = -U„( Z ) 

hn 



x n (z) = z n (l + zY 



n < 



Now it is easy to see the reason why Vn (z) does not satisfy (B.ll) in the generic situation. This is since the series 

-1 



(B.13) 
(B.14) 

(B.15) 



is divergent at small z though formally it is e qual to —1 and this formal answer miraculously coi ncide s with t he exact 

(2) \ k 1, V 

expression which can be obtained for V n (z) from (B.3) after summation from k = to k = n and using (B.13), (B.14). This is 

a typical example of how the formal analysis leads to a wrong conclusion and the analytical consideration is of great importance. 
Instead of doing all this, one can prove in the same way that V^ 2) (z) = w n (z) (see (k.22h, (U.iol)) satisfies 



C nk T k {z) = z - 1 V n (z) 



where C nk is given by (4.18). It can be done by representing (3.71 in the form 



(B.16) 



(B.17) 



( 2) 

and performing the summation from k = n to k = +00 . The only thing we need to prove is that V k (z) satisfies the 
boundary condition 



Ck(z 



)=±(^z- 1 V k+1 (z)-V l 



; 



+00 



(B.18) 



in the region of small 2 . This is in complete analogy with the proof of (B.4). 



The rest of equations in (1.13), (4.15) can be derived from (4.2J) using the same methods. 
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